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Interplay of kinematic and magnetic forcing in a model 
of a conducting fluid with randomly driven magnetohydrody- 
namic equations has been studied in space dimensions d > 2 
by means of the renormalization group. A perturbative ex- 
pansion scheme, parameters of which are the deviation of the 
spatial dimension from two and the deviation of the exponent 
of the powerlike correlation function of random forcing from 
its critical value, has been used in one-loop approximation. 
Additional divergences have been taken into account which 
arise at two dimensions and have been inconsistently treated 
in earlier investigations of the model. It is shown that in 
spite of the additional divergences the kinetic fixed point as- 
sociated with the Kolmogorov scaling regime remains stable 
for all space dimensions d > 2 for rapidly enough falling off 
correlations of the magnetic forcing. A scaling regime driven 
by thermal fluctuations of the velocity field has been identified 
and analyzed. The absence of a scaling regime near two di- 
mensions driven by the fluctuations of the magnetic field has 
been confirmed. A new renormalization scheme has been put 
forward and numerically investigated to interpolate between 
the e expansion and the double expansion. 

47.27.Gs,52.30.-q,11.10.Hi 



I. INTRODUCTION 

During the past two decades asymptotic analysis of 
stochastic transport equations [Navier-Stokes equation, 
magnetohydrodynamic (MHD) equations, advection- 
diffusion equation and the like] has attracted increasing 
attention. Various forms of the renormalization group 
(RG) have proved to be particularly useful in this in- 
vestigation, and a great deal of work has been carried 
out in the RG analysis of the stochastic Navier-Stokes 
equation and the problem of apassive scalar (turbulent 
diffusion or heat conduction) Somewhat less effort 

has been devoted to the asymtotic analysis of stochas- 
tic magnetohydrodynamics since the pioneering work of 
Fournier et al || and Adzhemyan et al [Q. In particular, 
in these papers the existence of two different anomalous 
scaling regimes (kinetic and magnetic) in three dimen- 
sions was established corresponding to two non-trivial 
infrared-stable fixed points of the renormalization group. 
It was also conjectured that in two dimensions the mag- 
netic scaling regime does not exist due to instability of 
the magnetic fixed point. However, in both papers there 
were flaws in the renormalization of the model in two 



dimensions (^J^]. Even more serious shortcomings are 
present in recent investigations of MHD turbulence HQ , 
in which a specifically two-dimensional setup has been 
applied with the use of the stream function and mag- 
netic potential. Therefore, results obtained for the two- 
dimensional case in these papers cannot be considered 
completely conclusive. 

In the present paper we have first carried out a field- 
theoretic RG analysis of the stochastically forced equa- 
tions of magnetohydrodynamics with the proper account 
of additional divergences which arise in two dimensions. 
This gives rise to a two-parameter expansion of scaling 
exponents and scaling functions the parameters of 
which are the deviation of the spatial dimension from two 
and the deviation of the exponent of the powerlike corre- 
lation function of random forcing from its critical value. 
In this double expansion the standard procedure of min- 
imal subtractions was used in the renormalization of the 
corresponding field-theoretic model. We have carried out 
a one-loop RG analysis of the large-scale asymptotic be- 
havior of the model and confirmed the basic conclusions 
of the previous analyses |||[ that near two dimensions 
a scaling regime driven by the velocity fluctuations may 
exist, but no magnetically driven scaling regime can oc- 
cur. We have also identified a scaling regime driven by 
thermal fluctuations || of the velocity field. 

Second, we have performed a renormalization of the 
model with a different choice of finite renormalization in 
order to find at which non-integer dimension the mag- 
netic fixed point ceases to be stable. This borderline 
dimension was found in Ref. Q with the use of the 
momentum-shell RG in a setup valid in a fixed space 
dimension d > 2. In the two-parameter expansion with 
the deviation of the exponent of the powerlike correlation 
function of random forcing from its critical value e and 
25 = d — 2 as expansion parameters this effect cannot 
be traced. Therefore, we have carried out a RG analysis 
according to the "principle of maximum divergences" in 
the sense that we have included in the renormalization 
all graphs relevant in two dimensions, and fixed the finite 
renormalization in a way which reproduces the results of 
a momentum-shell renormalization of Wilson J^] at one- 
loop order. This procedure gives rise to RG functions 
such that in the limit of small S, e they reproduce the 
results of the two-parameter expansion, and in the limit 
of small e (but finite 5) they yield the results of the usual 
e expansion . 

We have also investigated the long-range asymptotic 
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behavior of the model in the framework of the latter 
scheme without any small parameter and found, in par- 
ticular, that in this case thermal fluctuations make the 
value of the borderline dimension of the magnetic scal- 
ing regime significantly lower (d c = 2.46) than in the e 
expansion ||) (d c — 2.85). 

This paper is organized as follows: Section II starts 
from the functional formulation of the solution of stochas- 
tic MHD. This is convenient for the analysis based on 
the standard field-theoretic RG approach, the details of 
which are described in Sec. III. The Kolmogorov con- 
stant for MHD is calculated in Sec. IV at the leading 
order of the two-parameter expansion. In Sec. V a RG 
analysis of the model with maximum divergences is car- 
ried out in arbitrary dimension, and the observed strong 
effect of thermal fluctuations is discussed. In Sec. VI the 
conclusions are presented. 



II. FIELD THEORY FOR STOCHASTIC 
MAGNETOHYDRODYNAMICS 

We consider the model of stochastically forced con- 
ducting fluid described by the system of magnetohy- 
drodynamic equations for the fluctuating velocity field 
v(t,x) = v(ie) of an incompressible conducting fluid and 
the magnetic induction B = yfpjlh (p is the density and 
\x the permeability of the fluid) in the form 

d t v + P [(v • V)v - (b • V)b] - *A)V 2 v = P , (2.1) 
d t b + (v • V)b - (b • V)v - zy woV 2 b = f b , (2.2) 

together with the incompressibility conditions 

V • v = 0, V • P = 0, V • P = 0. (2.3) 



In (2.1) P is the transverse projection operator, uq the 
(unrenormalized) kinematic viscosity and 1/uq the un- 
renormalized magnetic Prandtl number. In statistical 
applications of the field-theoretic RG the unrenormalized 
(bare) paramaters are the physical ones. 

The statistics of v and b are co mple tel y d etermined 
by the non-linear equations (2.1), ( |2.2| ), ( |2.3D , and the 
probability distribution of the external large-scale ran- 
dom forces P, P. It is customary j3|Q| to consider ran- 
dom forces P and f b having a zero-mean Gaussian dis- 
tribution with correlation functions of the form 



Ann 0*0 = S(t) 



{2n) d 



e lk *D(k) S. 



k 2 



(2.4) 



in which the time correlations of the fields have the char- 
acter of white noise, while the spatial correlations are 
controlled by t he s calar function D(k). Transversality 
of the matrix (2.4) is a consequence of the equations 
V • v = V • b = 0. 



To analyze renor maliz at ion ne ar t wo dimensions we 
consider the model (2.1), (|2.2| ), (|2.3| ) supplemented by 
the forcing statistics 



[ fm( X l)fn( X 2) > = u ^0 D mn {X\ ~ X 2 ] {1, 9vW, 9v2o\) 

[f b m (xi)f b n (x 2 )) 



u v l D mn (Xl - X 2 \ {a, 5610, ffb2o}) , 



= 0, 



(2.5) 



where 



d d k 



Dmn (x; {A, B, C}) = 6{t) 
P mn (k)e lk - X [Bk 2 - 2S - 2Ae + Ck 2 ] . (2.6) 



All the dimensional constants g v w, <?bio, gv20 and gbw, 
in ( |2.5| ) control the amount of randomly injected energy. 
The choice of the values of the parameters e and S de- 
termines the powerlike falloff of the long-range forcing 
correlations and the space dimension of the system un- 
der consideration. 

We choose uncorrelated kinematic and magnetic driv- 
ing [(Pf b ) = 0], because we are considering arbitrary 
space dimension d > 2 and it is not possible to define a 
nonvanishing correlation function of a vector field and a 
pseudovector field in this case. This can be done sepa- 
rately for integer dimensions of space, but, contrary to 
claims of some authors ||,[l(| , is no obstacle for applica- 
tion of the RG H . 

The structure of the matrix D mn in (fOj) reflects a 
more detailed intrinsic statistical definition of forcing, 
whose consequences are deeply discussed in Refs. P,pT|- 
Technically, it is necessary to accompany the long-range 
corr elations [corresponding to the term / c 2 -2«5-2Ae m 
( |2.6| )] by local correlations [described by the analytic in 
k 2 term in the correlation function (2.6)] in order to con- 
struct a consistent renormalizationprocedure for the cor- 
responding field-theoretic model [p| JT^Jl^[ . This feature 
has been overlooked in the previous analyses of the prob- 
lem IHH] . The prefactors uq Vq and Uq v$ in (21) have 
been extracted fo r th e convenience of calculations. 

The definition (2.6) includes two principal - low- and 
high-wave-number scale - kinetic forcings separated by 
transition region in the vicinity of the characteristic wave 
number of order 0([S/C] 1 ^ 2l5+2Ae) ). The forcing con- 
tribution with local correlations gives a phenomenologi- 
cal description of small-scale thermal fluctuations of the 
magnetic induction and the velocity field [||. 

The long-range part s of the translational invariant cor- 
relation functions (2.6) become scale invariant at the val- 
ues e = 2, a = 1. For the exponent e the value e = 2 
is physically most reasonable, since it represents the as- 
sumption that random forces in the Navier-Stokes equa- 
tion (2.1) act at very large scales, which substitutes for 
the effect of boundary conditions. 

We are working in an arbitrary dimension, but the 
renormalization will be carried out within the two- 
dimensional model. In two-dimensional magnetohydro- 
dynamic turbulence, in contrast to fluid turbulence, there 
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are direct energy cascades in both two and three dimen- 
sions. Therefore it is natural to expect that the scaling 
behavior is rather similar in both cases, and we apply the 
same forcing spectrum in all space dimensions d > 2. 
We use the correlation functions 



( v jl ( x l) v 32 (X2)v h (x 3 ) ■ ■ ■ V jN (X N ) ) 



(2.7) 



where 1 < j r < d, r ~ 1, 2, • • • N as measurable quanti- 
ties for the description of turbulence statistics. We have 
applied the RG method to the calculation of asymptotic 
properties of the correlation functions in the way initi- 
ated in Ref. |Q. This approach is based on a formal 
mapping from the stochastic model (2.4) to a quantum- 
field model |l5|,|l6| with a De-Dominicis-Janssen action 
S{v, v', b, b'}, which is a functional of the physical fields 
v, b and independent solcnoid al au xiliary fields v', b'. 
Thus, the correlation functions ( |2.7| ) can be expressed as 
functional averages with the "weight functional" W = 
exp S. Th e sys te m of the stochastic MHD equations 
(|T]), (U), Q, (|U), and gives rise to the fol- 

lowing Dc-Dominicis-Janssen action: 

S = ~ Jdxi J dx 2 

x uqi/q v' m (x{)D mn {x\ - x 2 ;{l,g vl0 ,g v20 })v' n (x 2 ) 

+ u 2 l is a 3 b' m {x 1 ) D mn (xi - x 2 ;{a,g bWl g b20 }) b' n (x 2 ) 

+ J Ax jv' • [-d t v - (v • V)v + ^ V 2 v + (b • V)b] 

+ b' • [-d t b + u v V 2 h + (b • V)v - (v • V)b] } . (2.8) 

The dimensional constants g v io, g b io, g V 20, and g b 20, 
which co ntro l the amount of randomly injected energy 
through (2J)), play the role of expansion parameters of 
the perturbation theory. 



III. TWO-PARAMETER EXPANSION OF THE 
MODEL 



The action ( |2.8| ) gives rise to four three-point interac- 
tion vertices defined by the standard rules 0, and the 
following set of propagators 



A™„(k, t) = A^(-k, -t) = 0(t) P mn (k) e 



A^(k, t) = A^„(-k, -t) = 6{t) P mn (k) , 
A™„(k,i) = i Uo ^ 2 P m „(k)e^ fc2 l t l 



-Vq k t 
— UQU^k 2 t 



x (g vW k 



-2e-2S 



gv2o ) , 



(3.1) 



x (g b iok~ 2at - 25 



5f)20 J 

in the time-wave-number representatio n. W ith due ac- 
count of Galilei invariance of the action (2.8), and careful 



analysis of the structure of the perturbation expansion it 
can be shown ji} that for any fixed space dimension d > 2 
only three one-particle irreducible (IPI) Green functions 
put) ^ Y bb and T v hb with superficial UV divergences are 
generated by the action. They give rise to counter terms 
of the form already present in the action, which thus 
is multiplicatively renormalizable by power counting for 
space dimensions d > 2. 

We would like to emphasize that the structure of renor- 
malization should always be analyzed separately and it 
is not at all obvious that the nonlinear terms are not 
renormalizcd in the solution of the stochastic MHD equa- 
tions. In fact, direct calculation shows that the Lorentz- 
force term is renormalized. There seems to be a cer- 
tain amount of confusion about this point in the recent 
literature. For instance, in Refs. |p|JlC|l the authors er- 
roneously neglect renormalization of nonlinear terms as 
high-order effect. The approach adopted in Ref. || for 
two-dimensional MHD turbulence was quite recently crit- 
icized by Kim and Young |Q, who, alas, in their field- 
theoretic treatment of the same problem ignore renor- 
malization of the Lorentz force without any justification. 
They also neglect renormalization of the forcing corre- 
lations by effectively considering renormalization of the 
model at d > 2, which does not seem to be appropriate in 
a setup in which the strictly two-dimensional quantities, 
the stream function and magnetic potential, are used for 
the description of the problem. 

The analysis of the autocorrelation functions of the 
velocity field and magnetic induction is essential near 
two dimensions, since in two dimensions additional di- 
vergences in the graphs of the IPI Green functions T v v 
and r b b occur. The point here is that the nonlocal 

term of the action (and the similar one with the auxiliary 
field b') 

J At Jd d Xl y"d d x 2 v'(x 1 ,<)-v'(x 2 ,t) 



d d k 

(2n) d 



£,2-2<5-2e e ik-( Xl -x 2 ) 



brought about by the force correlation functions (2.5) 
is not renormalized since the divergences produced by 
the loop integrals of the graphs are always local in space 
and time |M. The simplest way to include the corre- 
sponding local counter terms v'V 2 v' and b'V 2 b' in the 
renormalization is to add corresponding local terms to the 
force correlation function at the outset in order to keep 
the model multiplicatively renormalizable, which is con- 
venient technically. This is why we have used the force 
correlation functions (2.5) and ( |2 . 6| ) with both long-r ange 
and short-range correlations. As a result, the action ( |2.8| ) 
is multiplicatively renormalizable and allows for a stan- 
dard RG asymptotic analysis jlTj. 

In the momentum-shell analysis of Fournier et al. jj| 
these divergences were taken into account only in the 
special case, when the force correlation function (|2.6| ) is 
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local, i.e. oc k 2 (formally this was fixed by the condi- 
tion 28 + 2Ae = 0, which excludes At from the param- 
eters of the model). In the field-theoretic treatment of 
Adzhemyan et al. the contribution of the additional 
divergences was prescribed to a renormalization of the 
non-analytic term cx /c 2 ~ 2t5 ~ 2e , although only analytic in 
k 2 terms are produced in the course of renormalization. 

The model is regularized using a combination of ana- 
lytic and dimensional regularization with the parameters 
t and 26 = d — 2. As a consequence, the UV divergences 
appear as poles in the following linear combinations of 
the regularizing parameters: t, S, 2t + d, and (a+ l)e + <5. 
The UV divergences can be removed by adding suitable 
counterterms to the bas ic action Sb obtained from the 
unrenormalized one (2^8 ) by the substitution of the renor- 
malized parameters for the bare ones: 3^10 — > /i 2e 3„i, 
g V 20 — > fi~ 2S g V 2, ffbio — > (t 2ae gbl, 9b20 — > M _2<5 5b2, Vq — > v, 
uo — > m, where /U is a scale setting parameter having the 
same canonical dimension as the wave number. 

To construct an analog of the usual e expansion with 
e and 5 as small parameters of the same order of mag- 
nitude, it is convenient to use the minimal-subtraction 
(MS) scheme for the renormalization. In this approach 
only singular parts of the Laurent series of the superfi- 
cially divergent 1PI Green functions are included in the 
renormalization constants, which give rise to the counter 
terms added to the basic action to make the Green func- 
tions of the resulting renormalized model UV finite. The 
counter terms for the b asic action corresponding to the 
unrenormalized action (2.8) are 



AS = J dx [v (Z x - 1) v'V 2 v + uv (Z 2 - 1) b'V 2 b 



+ -(l-ZM„2A" M v'V 2 v' 

+ 1(1- Z 5 )u 2 v 3 g b2f i- 2S b'V 2 b' 
+ (Z 3 -l)v'(b-V)b] , 



(3.2) 



where the renormalization constants Z\, Z 2 , Z&, Z§ 
renormalizing the unrenormalized (bare) parameters 
eo = {3^10, 9v20, 9bw, gb2Q,Ua,Vo} and the constant Z 3 
renormalizing the fields b, and b', arc chosen to can- 
cel the UV divergences appearing in the Green functions 
constructed using the basic action. Due to the Galilean 
invariance of the action the fields v', and v are not renor- 
malized. 

In a multiplicatively reno rma lizable model, such 
as (pTsJ) , the counter terms ( |3.2| ) can be chosen in a 
form containing a finite number of terms of the same al- 
gebraic structure as the terms of the original action ( 2 



Thus, all UV divergences of the graphs of the perturba- 
tion expansion may be eliminated by a redefinition of the 
parameters of the original model. 

Renormalized Green functions are expressed in terms 
of the renormalized parameters 



'Jvl 



gvio A* 2e Z 2 Z 2 , g V 2 



9v20 (J 26 ZfZ 2 Z 4 1 , 



U0 Z-y , u = ixq Z2 ^ Z\ , 

— 2 ae r? ry 2 ry — 1 

Zj\Zj 2 Zj^ 



V = 

'61 = 3610 A* 



362 — 3620 M 



(3.3) 

2<5 ry ryl ry — 1 ry — 1 
Zj\Zj 2 /j z A 5 , 



which are the parameters of the renormalized action 
Sr = Sb + AS connected with the unrenormalized action 
( |2.8| ) by the relation of multiplicative renormalization 

S R {v, b, v', b', e} = S {v, bzjy, h'Zp, e } , 

where e is a shorthand for all the renormalized param- 
eters 3„i, 3„2, 36i, 362, u, and v. Calculation of the cor- 
relation and response functions of the velocity and mag- 
netic fields with the use of the renormalized action yields 
renormalized Green functions without UV divergences. 

Independence of the unrenormalized Green functions 
of the scale-setting parameter [i may be expressed in the 
form of differential RG equations for the renormalized 
1PI Green functions. To keep notation simple, we quote 
these equations only for the renormalized pair correlation 
functions of the velocity field and the magnetic induction. 
We define the Fourier transforms as 



W v R v mn (h 



*2,k;g) 



d d Xl 



(2ir) d 

x( Um (x 1 ,i 1 )«„(x 2 ,t 2 ))e 4k -( Xl - X2 ), 



Wi 



66 

R mn 



*2,k;3) = / j^ 1 



(2n) d 

x(fe m (x 1 ,i 1 )fe„(x 2 ,< 2 ))e 4k -( Xl - X2 ). 
The basic RG equations for the correlation functions are 



ft 5 



7i v 



dv 



W. 



R mn 



o. 



where (3 g d g is a shorthand for 



W. 



66 

R mn 



= o, 



(3.4) 



z Pgvl-^ \-0gv2 7; \~Pgbl 



dg v i 



' dg v2 



dg, 



61 



~Pgb2 



dg, 



-0v 



1,2 



du 



The coefficient functions of Eqs. (3.4) (3 g and 71 are ex- 
pressed in terms of logarithmic derivatives of the renor- 
malization constants. We use the definitions 



d\nZi 



Of 1 



a d 9 



(3.5) 



where 3 = {3^1, 3^2, 361, 362, u}, and the subscript "0" 
refers to partial derivatives taken at fixed values of the 
bare parameters eo- It should be noted that here the 
functions f3 g and 71 are functions of the parameters 3 
only. 

Expressing the correlation functions through dimen- 
sionless scalar functions R v and Rb as 
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W v R v mn (t,k;g) = v 2 k- 25 P mn (k)R v (T,s;g) , 
W b R b mn (t,k;g) = V 2 k- 2S P mn {k)R b (T,s-g), 

where s = k/n, s € [0, 1] is the dimensionless wave num- 
ber, and t = tvk 2 the dimensionless time, and solving 
Eqs. ( |3.4[ ) by the method of characteristics, we obtain 
the correlation functions in the form 

W% mn {t,k;g) = P mn {k)V 2 k~ 25 R v (tk 2 V,l;g) , 

W bb mn {t,k;g) = J>^^ (3.6) 
x P mn (k)V 2 k- 2S R b (tk 2 V,l;g) , 

where g is the solution of the Gell-Mann-Low equations: 

dg(s) 



din. 



09 W)\ 



(3.7) 



and v is the running coefficient of viscosity 
77 = ve ~ I" Ax ^^))/^ ^ 

The scale-invariant asymptotic behaviour of the correla- 
tion functions stems from the existence of a stable fixed 
point of the RG transformation @ g = determined by 
the Gell-Mann-Lo w equ ations (3.7). 



The definitions (3.5) and the relations (3^) yield (3 
functions of the form 



0gvi = 9vi (-2e + 271 + 72) , 

0gv2 = 9v2 (2(5 + 27i + 72 - 74) , 

g bl = .961 (-2 ae + 71 + 272 - 73) , 

0gb2 = 9b2 (2(5 + 71 + 272 - 73 - 75) , 
0u = U (71 - 72) . 

At one-loop accuracy the 7 functions are 

7i = -r~[u (9vi + 9v2) + 9bi + 9b2] , 

OZ7T 



(3.8) 



72 = 

73 = 

74 = 



1 9vl + 9v2 - gbi - 9b2 



1 



8tt 
1 

1 u (g vl + g v2 ) 2 + {g b i + flto) 2 



{9bi + 9b2 - 9vi ~ g V 2) , 



(3.9) 



32tt 



gv2 



There are no UV divergences in the 1PI Green function 
r b b in the one-loop approximation, therefore 



75 = 0, 



1. 



(3.10) 



which is a specific property of the two-dimensional MHD. 

Large-scale asymptotic behavior is governed by in- 
frared stable fixed points of Eqs. (3.7), determined by 



the system of equations g {g*) = 0, and the conditions 
g — > g* , when s — > 0. For g(s) close to g* we obtain a 
system of linearized equations 



Is- CI 

as 



(9-9*) 



where I is the 5x5 unit matrix and the matrix f2 = 
(d[3 g /dg)\ g > . Solutions of this system behave like g = 
g* + 0(s j ), when s — » 0. The exponents \j, j — 
1,2,3,4,5 are the eigenvalues of the matrix CI. In the 
vicinity of the fixed point all the trajectories g(s) ap- 
proach the fixed point, if the matrix CI is positive definite 
[i.e. Re(Aj) > 0]. 

Apart from the Gaussian fixed point g*-^ = g* 2 = g b i = 
■9b2 = ^ with no fluctuation effects on the large-scale 
asymptotics, which is IR stable for(5>0,e<0,a>0, 
there are two nontrivial IR stable fixed points of the RG 
with nonnegative g* 1: g* 2 , g^ v g^ 2 , and u*. 

The thermal fixed point is generated by short-range 
correlations of the random force with 



0. 



g* v2 = -4tt(1 + V17)S, 

gti = °> s&2 = °> 



and the inverse magnetic Prandtl number 



Vvf-i 



1.562. 



(3.11) 



(3.12) 



Physically, the asymptotic behavior described by this 
fixed point is brought about by thermal fluctuations of 
the velocity field [p|. Th e reg ion of stability of the ther- 
mal fixed point fl3.1l| ), ( |3.12| ) is 2e + 3<5 < 0, S < in 
the 5, e plane. For the magnetic forcing-decay parame- 
ter a the stability region is determined by the inequality 
8ae + (13 + Vl7)S < 0. 

The kinetic fixed point [g| generated by the forced fluc- 
tuations of the velocity field is given by the universal 
inverse magnetic Prandtl number (3.12), the parameters 



.9,1 



128 7T e(2e + 3<5) 



9(VT7-1) e + 5 



g V 2 



1287T 



(3.13) 



9(v / 17- 1) S + e' 
and zero couplings of the magnetic forcing 

9bi = .962 = , 

and it may be associated with turbulent advec tion o f the 
magnetic field. The values of g* x and g^offi ( 3.13 ) cor- 
respond to those find previously in Ref. |5|. The region 
of stability of the kinetic fixed point in the 5, e plane 
is e > 0, 2e + 3(5 > 0. The stability of this fixed point 
also requires that the parameter a < (13 + VT7)/12 w 
1.427 independent of the ratio 5/e. In spite of the ab- 
sence of renormalization of the forcing correlation, the 
momentum-shell approach Q yields the same condition. 

The system of equations d3.8| ), ( |3.9| ), and (3.10) for 
the fixed points in this multi-charge problem is rather 
complicated and thus has several (in general complex- 
number) solutions, which we do not quote explicitly here, 
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because they are not physically relevant: apart from 
the fixed points listed above there are eight IR unsta- 
ble real-number fixed points in the physical region (all 
g > 0) of the parameter space, and several unphysical 
ones. Among the unstable fixed points are, in particu- 
lar, all the possible candidates to magnetic fixed points, 
i.e. fixed points with a nonvanishing magnetic coupling 
constant. Therefore, the conclusion made in Refs. [g|JJ] 
(although on inconsistent grounds) that the RG does 
not predict any magnetically driven scaling regime at 
and near two dimensions, is confirmed in the double- 
expansion approach. 

It should be noted, however, that the vanishing of the 
function 75 renders the linear combinations of 7 functions 
in the functions (3 g bi and /3 g b2 e qual: they both contain 
only 71 + 272 — 73 [see Eq (3.8)]. This has the impor- 
tant consequence that there is no fixed point with both 
gbi and gbi nonvanishing: at least one of them must be 
zero. This, of course, severely reduces the set of possible 
magnetic fixed points at the outset. 

It would be interesting, however, to follow the 
crossover from this regime to the scaling regime gov- 
erned by the competition of the stable kinetic and mag- 
netic fixed points, which exists in three dimensions. In 
the double-expansion approach the space dimension is as- 
sumed to be close to two, therefore the results obtained 
above are not applicable to this end. In the usual e expan- 
sion the leading-order value of the borderline dimension 
between the two regimes d c = (3 + Vfi49) /10 « 2.848 [§ . 

All the renormalization constants and the RG func- 
tions quoted above may be calculated also at finite S. The 
resulting system of fixed-point equations allows for a solu- 
tion in a form of an e expansion (with finite 5) and yields 
the same result as the usual e expansion at the leading 
order. However, this approach is not self-consistent in 
the sense that the field theory is not renormalizable at 
finite 5, but only in the form of a simultaneous expansion 
in the coupling constants and <5 Jl7[ . Therefore, in order 
to construct an interpolation procedure between the e, <5 
expansion and the usual e expansion something has to be 
done with the UV divergences at d > 2 introduced by the 
local terms in the force correlation functions. This issue 
will be dealt with in Sec. V. 



IV. KOLMOGOROV CONSTANT FOR 
STOCHASTIC MAGNETOHYDRO DYNAMICS 
NEAR TWO DIMENSIONS 



we infer an expression for the energy spectrum in terms 
of the scaling functions as 

(d-l)k^ 3 ( g vl0 u Q ^\ 2/3 



E{k) (47r)<V2r(d/2) V g vl u 

R v (0, l;g) + Z 3 (g)ei: dx ^ x))/x R b (0, l;g) 



where for the running coefficient of viscosity v the ex- 
pression 



3\ 1 / 3 



(4.2) 



has been used. The relation (4.2) is a consequence of the 
connection between the functions (3 U , 71 and 72 [ p"8[ . 
At the kinetic fixed point the spectrum has the form 



[ ' { g* vl u* J (An)^T(d/2) 
R v (0, !;<?*) + Z 3 {g)s^ 3 ^R b (0, !;<?*) . 



(4.3) 



However, at leading order i?&(0, l;g) = \{gbi + gbz), and 
since at the kinetic fixed point <j£i = 9b2 = an d 73 (<7*) > 
0, only the scaling function R v survives here. 

The kinetic and magnetic energy injection rates e v and 
Eb may be expressed as 



Ey 
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1 
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d rf k 
(2^p ' 
d d k 
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;r(k)-p(-k)), 



;f b (k).f b (-k)>. 



(4.4) 



For our choice of correlation functions after the introduc- 



tion of simple sharp cutoffs, Eq. (4.4) yields the relation 
between the unrcnormalized values of the coupling con- 
stants and the energy injection rates in the form 
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(4.5) 



The energy spectrum E(k) in magnetohydrodynamics 
is given by the sum of equal-time pair correlation func- 
tions of the velocity field and magnetic induction: 

W™(t,x;i,x) + tC(i,x;i,x) = 2 / dk E(k) . (4.1) 

Jo 



Note that the energy spectrum is defined through un- 
renormalized correlation functions. From (|3.6|) and (4.1) 



where kj is the wave number corresponding to the inte- 
gral scale, and kd, k' d are the characteristic wave numbers 
of viscous and resistive dissipation, respectively. 

In the stationary state modeling developed isotropic 
turbulence, the energy injection is assumed to take place 
at large scales. Therefore, we put the parameters g V 2o 
and <?fc20i which correspond to small-scale injection of en- 
ergy, equal to zero. It should be borne in mind that the 
corresponding running coupling constants are created in 



G 



the course of renormalization regardless of the unrenor- 
malized values of these parameters. The present pertur- 
bative calculation yields only the leading order in the e, 
8 expansion of the amplitude coefficients in the scaling 
form of the correlation functions. The refore the coupling 
constants should be solved from (4.5) as functions of e v 
and Eb also only at leading order of e, 8 expansion. Thus, 
we arrive at the relations 



_ UqVq g vW 4 _ 2e 
~ 167T d ' 

167T {kd) 



(4.6) 



valid for large Reynolds and magnetic Reynolds numbers, 
when kd/kj ~ Re 3/4 > 1 and kUkj ~ Rm 3/4 > 1. 



Substituting the relations (4.6) in Eq. (4.3) we arrive 
at the spectrum in the form 



(u*) 1/3 (g*i +g* v2 ) 2/3,a-4e/3i.4(6-2)/3 
(2^)1/3(3^)2/3 ^ 

= C fe ^/3 fc l-4e/3 fc 4(e-2)/3 



at the leading order of the 8, e expansion. The value of 
the Kolmogorov constant Ck inferred from here 



Cfc = 



2 ■ 12 1 ' 3 e 1 ' 3 (e + S) 2 ' 3 
(2e + 3<5) 2 / 3 



coincides with that obtained in the case of turbulent ad- 
vection of a passive scalar Jl^ ], which is, of course, not 
surprising, since the magnetic field is passively advected 
by the velocity field in the scaling regime governed by 
the kinetic fixed point. 



V. RENORMALIZATION WITH MAXIMUM 
DIVERGENCES ABOVE TWO DIMENSIONS 

We want to maintain the model UV finite for 28 = 
d—2 > and simultaneously keep track of the effect of the 
additional divergences near two dimensions. To this end 
we introduce an additi ona l UV cutoff in all propagators, 
i.e. instead of the set (3.1) we use the propagators 



A™ n A (k, t) = 9{t) 0(A - k) P mn (k) e-»° k 1 , 



A b „ b ^(k,t)=6(t)e(A-k)P mn (k)e 



—uquqU t 



C(k, t) = - 6(A - k) u a 14 P mn (k) e 



-v k |t| 



(gvio k 



-2e-25 



9v20 , 



A^(k, t) = - 6(A - k) u a 14 P mn {k) e 



-UQi^ok \t\ 



(gbw k 



-2ae-2S 
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where A is the cutoff wave number. This change ob- 
viously does not affect the large-scale properties of the 
model. We would like to emphasize that the additional 



cutoff must be introduced uniformly in all lines in or- 
der to maintain the model multiplicatively renormaliz- 
able. An attempt to introduce the cutoff, say, in the 
local part of the correlation functions only by the sub- 
stitution k 2 — * 6 (A — k)k 2 would fail to renormalize the 
model multiplicatively, because loop contributions to the 
complete (dressed) correlation function would not repro- 
duce such a structure in the wave-vector space. 

In contrast with particle field theories we will keep the 
cutoff parameter A fixed, although large compared with 
the physically relevant wave-number scale. This intro- 
duces an explicit dependence on A in the coefficient func- 
tions of the RG, which has to be analyzed separately in 
the large-scale limit in the coordinate space. The setup 
is thus very similar to that of Polchinski J|(| . 

The RG equations maintain the previous form (|3.4| ), 
but the coefficient functions become, in general, functions 
of the parameters n and A through the dimensionless 
ratio /z/A. Solving the RG equation by the method of 
characteristics we obtain the solution 

W^ mn (t,k, A;g) = P mn {k)V 2 k- 2S R v (tk 2 V, 1, ^-;g 



W. 



bb 

Ft mn 



(t,sk,A;g) = ef" dxi3ix)/x 



\ />, k)V-k---'fh[tk 2 V,l,!j : ;g) , 



where g is now the solution of the Gell-Mann-Low equa- 
tions: 



rflns 



sO)) 



(AS 

T 



with the [3 functions explicitly depending on s, the di- 
mensionless wave number. 

Above two dimensions an UV renormalization of the 
model would require and infinite number of counter terms 
and in this sense it is not renormalizable in the limit A — > 
oo. To avoid dealing with these formal complications, we 
keep the additional cutoff A fixed (although large), and 
choose the renormalization procedure according to the 
principle of maximum divergences Jis[| : the same terms 
of the action are renormalized as in the two-dimensional 



case in the previous section (3.2), but the renormaliza- 
tion constants may have an explicit dependence on the 
scale-setting parameter through the ratio (J./ A. The two 
counter terms / dec [| (1 — Z i )uu 3 g V 2 H~ 25 v'V 2 v' + 1 (1 — 
Z 5 ) u 2 v 3 gb2 fJ-~ 25 b'V 2 b' ] are superfluous in the sense 
that in the limit /x/A — > the contribution to the Green 
functions of the graphs containing the coupling constants 
g V 2 and gb2 remains finite provided 28 — d — 2 is fixed 
and finite and the other counter terms are properly cho- 
sen. This is guaranteed by Polchinski's theorem |20|. We 
retain these counter terms in order to have a possibility 
to pass to the limit 8 — » in the RG equations. 

The presence of the explicit cutoff implies some tech- 
nical difficulties in the calculation of the renormaliza- 
tion constants in the traditional field-theoretic approach, 
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which arise because we are dealing with massless vec- 
tor fields. It turns out that the coefficient functions of 
the RG equation are simplest in a renormalization pro- 
cedure, which is similar to the momentum-shell renor- 
malization of Wilson H . If we were calculating over the 
whole wave-vector space without an explicit UV cutoff, 
there would not be much difference between the effort re- 
quired in both approaches. The presence of the UV cutoff 
makes calculations with nonvanishing external wave vec- 
tors rather tedious. 

Let us remind that the choice of a renormalization pro- 
cedure basically is the choice of the rule according to 
which the counter-term contributions are extracted from 
the perturbation expansion of the Green functions of the 
model. The usual field-theoretic prescription goes as fol- 
lows |2^|: Consider a 1PI graph 7, let R(j) be the renor- 
malised value of the graph, and let R("f) the value of the 
graph with subtracted counter terms of all the subgraphs, 
then 



i?( 7 ) = R(l) - TR(>y) 



(5.1) 



where the operator T extracts the counter-term contri- 
bution from i?(7). Usually R(^) = 7 on one-loop graphs, 
and the renormalization scheme is specified by the ac- 
tion of R(j) on multiloop graphs together with the def- 
inition of the operator T. The counter ter ms m ay then 
be constructed recursively with the use of (5.1) and the 
definitions of T and i?( 7 ). There is rather large freedom 
in the choice of the counter-term operator, but to ar- 
rive at Green functions finite in the limit A — > 00 in two 
dimensions - which we want to have a connection with 
the double expansion - the operation R must be chosen 
properly. 

Here, we have used a renormalization procedure, in 
which the operation i?( 7 ) is standard and the sub- 
traction operator T is defined as follows: let F 7 (w,k, A) 
be the function of external frequencies and wave-vectors 
(which also depends on the cutoff parameter A) corre- 
sponding to the expression i?( 7 ) (this is not a 1PI graph, 
in general). The subtraction operator T extracts the 
same set of terms of the Maclaurin-expansion in the ex- 
ternal wave-vectors, which generate the counter terms 
(3.2), from the difference F^(u, k, A) — F 7 (w, k, fi). These 
coefficients of the Maclaurin-expansion are calculated 
at vanishing external frequencies and wave-vectors. It 
should be noted that the coefficients of this Maclaurin 
expansion of the function -F 7 (w, k, A) itself may not exist 
in the limit u> — ► in this "massless" model, but the dif- 
ference .F 7 (w,k, A) — Fj(u>, k, ff) allows for a Maclaurin 
expansion finite in the limit u) — > to the order required 
for the renormalization. The counter-term operator T 
and the combinatorics of the renormalization procedure 
for higher-order graphs may then be constructed in the 
standard fashion. Although this is actually not needed in 
the present one-loop calculation, the very possibility of 
this extension is necessary to guarantee that the renor- 
malization renders the model finite in the limit A — > 00 
in two dimensions. 



Effectively, at one-loop order this prescription reduces 
the region of integration to the momentum shell [i < 
k < A, which leads to the same calculation as in the 
momentum- shell renormalization. In higher orders, how- 
ever, our renormalization scheme does not coincide with 
the momentum-shell renormalization. The point of the 
present renormalization procedure is that without some 
sort of IR cutoff the subtraction at zero momenta and fre- 
quencies is, in general, not possible in a massless model, 
whereas a subtraction at vanishing frequencies and ex- 
ternal momenta of the order of fi leads to much more 
complicated calculations due the heavy index structure. 

At one-loop accuracy in this scheme the 7 functions 
are 



7i 
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2B 



(d 2 -d- 2e) ug v i + (d 2 + d - 4 + 2ae) g bl 

+ (d 2 - 2) (u g v2 +g b2 ) , 
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l+u)B 



(d-l)(d + 2) (ffvi+flta) 



{d + 2) (d-3) (g b i+g b2 ) 
2 
B 

d 2 -2 



[gbi + gt2 - g v i - 3^2] 



(5.2) 



2g v2 B 
2(d-2)(d+2) 
g b2 (l + u)B 



u (9vi + 9v2) 2 + (gbi + gb2) 2 



(gbi +gb2) (g v i + 5^2) 



where B = d(d + 2)T(d/2)(An) d / 2 . These expressions 
reveal an additional advantage of this renormalization 
scheme: at one-loop order there is no explicit dependence 
on /i/ A in the coefficient functions of the RG. At one-loop 
level a direct comparison with the expressions obtained 
in the framework of the Wilson RG is also possible: the 
dependence on g v \, g b \, and u of the f3 functions (3 gv \, 
(3 gb \ and (3 U corresponding to (|]^) coincides with that of 
their counterparts of Ref. || up to notati on. 

The set of (3 functions generated by (5.2) allows for 
a fixed-point solution in the form of an e expansion. 
Little reflection shows that the fixed-point equations in 
this case have a self-consistent solution with the fol- 
lowing leading-order behavior: u = 0(1), g v \ = O(e), 
.961 = 0(e), g v2 = 0(e 2 ) and g b2 = 0(e 2 ). The ac- 
tual fixed-point values of g v ±, g b \ and u in the e expan- 
sion as well as the stability regions with respect to e are 
determined by the same set of equations as in the ear- 
lier momentum-shell |Q and field-theoretic |Q calculation 
above two dimensions. The stability condition with re- 
spect to the dimension of space of these fixed points is, 
as expected, d > 2. 

It should be noted that the function 75 is finite in the 
set (^^), whereas in the double-expansion approach it 
was equal to zero [Eq. ( [3.10 )]. This means that magnetic 
fixed points with both g b i and g b2 may exist. In fact, 
there is one such fixed point stable at high dimensions 
of space which gives rise to a magnetically driven scaling 
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regime. This fixed point may be found in the e expansion, 
and we have also investigated its stability numerically. 
Technically speaking, the appearance of a magnetic fixed 
point with both magnetic couplings nonvanishing would 
be a completely expected thing to happen in the two-loop 
approximation, since we have not found any symmetry 
reasons or the like to prevent the renormalization of the 
magnetic forcing at higher orders. Thus, to investigate 
this effect consistently in the e expansion would require 
a full twodoop renormalization of the model, which is 
beyond the scope of the present analysis. 

On the other hand, in the case 6 ~ e — * the set of 



coeffi cient functions (5.2) yields the coefficient functions 
(3.9) of the double expansion of the previous section. 
Therefore, we think that it is not t otal ly unreasonable 
to use the set of coefficient functions (5.2) for an analysis 
of the RG fixed points for all dimensions d > 2. 

We have inves tiga ted the sta bility of the kinetic fixed 
point given by (3.8) and (5.2) numerically for a finite 
range of e with results depicted in Figs. [I] and |^. We 
think that this calculation, although it is a somewhat 
uncontrollable approximation, exhibits the effect of the 
thermal (short-range) fluctuations of the fields quali- 
tatively correctly. The stability of the kinetic scaling 
regime is strongly affected by the behavior of magnetic 
fluctuations: from Fig.[t] it is seen that the steeper falloff 
of the correlations of the magnetic forcing in the wave- 
vector space compared with that of the kinetic forcing 
the lower is the space dimension, above which the kinetic 
fixed point is stable. In particular, when the parame- 
ter a > 1.427 the kinetic fixed point ceases to be stable 
even in two dimensions. In three dimensions the kinetic 
scaling regime is stable against magnetic forcing, when 
a < 1.15. 

The monotonic growth of the kinetic-fixed-point value 
of the inverse magnetic Prandtl number u* as a function 
of the kinetic forcing-decay parameter in a fixed space 
dimension is depicted in Fig. |[ The plot shows also that 
u* is a monotonically decreasing function of the space 
dimension at fixed e. The lowest-lying curve corresponds 
to the leading order of the e expansion jjj 



8(d + 2) 



We are particularly interested in the stability of the mag- 
netic fixed point, and have carried out extensive numer- 
ical calculations of the stability of this fixed point as a 
function of e and the space dimension d. The results are 
plotted in Figs. || and |. 

In Fig. H the magnetic forcing-decay parameter a < 1 
(i.e. the kinetic-forcing correlations fall off steeper in the 
wave-vector space than those the magnetic forcing) and 
it is seen that for very small a a slowly enough decay- 
ing kinetic forcing renders the magnetic scaling regime 
unstable. In particular, this treshold is very small in 
three dimensions. With the growth of a a strip of stabil- 
ity of the magnetic fixed point in the e, d plane appears 



such that the magnetic regime remains stable in three 
dimensions for all allowed kinetic forcing patterns. It is 
also seen that the magnetic fixed point is persistently 
unstable at d < 2.46 for all e. This borderline dimen- 
sion should be compared with that given by the e expan- 
sion d c = 2.85. From the solution it can be seen that 
this significant discrepancy is due to the appearance of 
a stable magnetic fixed point completely different from 
that found in the e expansion: in the latter the magnetic 
fixed point is given by g*± = g* 2 = g£ 2 = u* = and 
flfo = U(d + 2)r(d/2)(47r) d / 2 ae/(d 2 - 3d - 32), whereas 
at the magnetic fixed point, whose stability is plotted in 
Figs. ^ and ^, only g* ± = u* = with nonvanishing fixed- 
point values of the other couplings. Thus, the lowering 
of the borderline dimension of stability of the magnetic 
scaling regime is a strong effect of the thermal fluctu- 
ations described by the short-range terms in the forcing 
correlation functions. Fig. ^ shows the lower boundary of 
the stability region of the magnetic fixed point for large 
values of a, when magnetic-forcing correlations fall off 
much faster that kinetic-forcing correlations in the wave- 
vector space. A remarkable feature of both plots is the 
insensitivity of the lower border of the stability strip to 
the value of magnetic forcing-decay parameter a. 



VI. CONCLUSIONS 

In conclusion, we have carried out a RG analysis 
of the large-scale asymptotic behavior of the solution 
of stochastically forced magnetohydrodynamic equations 
for all space dimensions d > 2. We have taken into ac- 
count the additional divergences appearing in two dimen- 
sions ignored or improperly treated in previous work. 
In a two-parameter expansion scheme we have found 
three infrared-stable fixed points in the physically rel- 
evant region of the parameter space spanned by the forc- 
ing parameters and the inverse magnetic Prandtl num- 
ber. Anomalous scaling behavior is brought about in 
the basins of attraction of the thermal and kinetic fixed 
points, the former of which is due to thermal fluctuations, 
and the latter due to long-range correlated random forc- 
ing of the Navier-Stokes equation. The thermal fixed 
point is related to the anomalous asymptotic behavior 
due to thermal fluctuations near two dimensions. With a 
proper choice of the force-correlation function the regime 
governed by the kinetic fixed point may be related to de- 
veloped isotropic turbulence in a conducting fluid. We 
have obtained the stability condition of the kinetic fixed 
point near two dimensions with respect to the magnetic 
forcing-decay parameter as a < 1.427, which coincides 
with that of Ref. [pf, but differs from the result of Ref. 

We have also put forward an interpolation scheme, 
which reproduces the earlier results in the case of an e 
expansion in any fixed space dimension d > 2, and the re- 
sults of the present paper in the two-parameter expansion 
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in e and 26 — d — 2 near two dimensions. Using this in- 
terpolation approach we have qualitatively analysed the 
dependence of the stability of the kinetic and magnetic 
scaling regimes on the forcing-decay parameters and the 
dimension of space and found that thermal fluctuations 
drastically lower the borderline dimension of stability of 
the magnetic scaling regime. 
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FIG. 1. The borderline dimension d c between the stability 
regions of the kinetic fixed point of the RG equations (3.5) 
and (^^) for magnetic forcing-decay parameter a near the 
double-expansion treshold a — 1.427. This plot reveals the 
strong dependence of the borderline dimension on the param- 
eter a. The shaded region on the right corresponds to values 
e > 2/ a, for whic h the forcing correlation function in the pow- 
erlike form ( |2.6[ ) leads to intractable IR divergences, and a 
corresponding IR cutoff (magnetic integral length scale) must 
be introduced. 
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FIG. 2. The fixed-point value of the inverse magnetic 
Prandtl number u* as a function of the space dimension d 
and the decay parameter e. The lowest curve corresponds to 
the leading order in the e expansion, which is not affected by 
thermal fluctuations. The shaded region in the upper part of 
the plot corresponds to values e > 2, for which an IR cut- 
off (kinetic integral length scale) must be introduced in the 
correlation function (2.C). 



FIG. 3. The borderline dimension d c between the stabil- 
ity regions of the magnetic fixed point of the RG equations 
( fo.8j ) and ( |5.2| ) for magnetic forcing-decay parameter a < 1. 
For sufficiently small values of a the magnetic fixed point is 
unstable for any finite value of e, but the region of stability 
grows with the growth of a so that for a > 0.658 the magnetic 
point becomes stable even in three dimensions for finite val- 
ues of e. The shading shows the re gion , where e > 2, in which 
the powerlike correlation function (2.6) cannot be consistently 
used. 
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FIG. 4. The borderline dimension d c between the stabil- 
ity regions of the magnetic fixed point of the RG equations 
( fo.g| ) and (^^) for large values of the magnetic forcing-decay 
parameter a > 1. The shaded area and half-planes with verti- 
cal dashed border lines show regio ns, w here e > 2/ a, in which 
the powerlike correlation function (2.6) cannot be consistently 
used. 
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